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Abstract: A Standard Model-like Higgs near 125 GeV in the MSSM requires multi-TeV 
stop masses, or a near-maximal contribution to its mass from stop mixing. We investigate 
the maximal mixing scenario, and in particular its prospects for being realized it in po- 
tentially reahstic GUT models. We work out constraints on the possible GUT-scale soft 
terms, which we compare with what can be obtained from some well-known mechanisms of 
SUSY breaking mediation. Finally, we analyze two promising scenarios in detail, namely 
gaugino mediation and gravity mediation with non- universal Higgs masses. 
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1 Introduction 



Recent results from ATLAS [1] and CMS [2] show intriguing hints for a Higgs boson with a 
mass around 124-126 GeV. While the reported excess of events is only at the level of about 
2-3 a per experiment, the consistency between the already excluded mass range, the excess 
in the remaining small window, and theoretical expectations provides a strong motivation 
to take these hints seriously and investigate their implications.^ 

If the cause of the observed excess were the lightest Higgs boson of the MSSM, it 
would be rather heavy, requiring large radiative corrections to its mass from the top-stop 
sector. In this case either the average stop mass must be large, Ms = ^ym^lTT^ > 3 TeV, 
or the stop mixing parameter \Xt\ must be around twice Ms [4]. The latter, known as the 
"maximal mixing scenario," is the subject of this paper. 

Several recent papers, including [5-13], have explored the implications of a heavy 
MSSM Higgs from a bottom-up perspective, prescribing the MSSM soft parameters at the 
TeV scale. ^ In this approach, the soft terms can simply be chosen by hand to yield maximal 
mixing. However, one should keep in mind one of the key motivations for low-energy 
supersymmetry: The supersymmetric Standard Model can naturally be extrapolated up to 
a very high fundamental scale. Indeed, gauge coupling unification in the MSSM points to 
the GUT scale Mqut ^ 10^^ GeV as the scale where it should be embedded into a more 
fundamental theory. It is therefore worthwhile to investigate if maximal mixing can result 
from some reasonable choice of GUT-scale parameters, what further relations between the 
GUT-scale soft terms this would imply, and which classes of models of high-scale SUSY 
breaking mediation can (or cannot) accommodate maximal mixing. Furthermore, it is 
clearly of interest how the GUT-scale conditions for maximal mixing affect the physical 
spectrum and the low-scale observables, such as Higgs cross-sections and decay rates. These 
subjects are addressed in the present paper. 

The implications of a 125 GeV Higgs for GUT-scale MSSM scenarios have been inves- 
tigated previously in [19-25].^ It was observed [19] (see also [31]) that, in the CMSSM and 
in NUHM models with sizeable mg, large l^ol ^ 2 mo is preferred to obtain a heavy Higgs. 
Our work goes beyond these studies by providing a thorough discussion of the prerequisites 
for maximal mixing in more general models, accompanied by a detailed numerical analysis. 

This work is organized as follows: In Section 2 we give a brief review of the maximal 
mixing scenario, and explain why it is non-trivial to obtain maximal mixing when running 
from the GUT-scale. Using semi-numerical solutions of the renormalization group equa- 
tions, we derive some necessary conditions for maximal mixing. In Section 3 we comment 
on the possibility of realizing maximal mixing in several well-established classes of models of 
SUSY breaking mediation, namely gaugino mediation, models with strongly-coupled near- 
conformal hidden sectors, radion mediation in 5D models, and gauge mediation. Section 4 
contains a detailed numerical analysis of a gaugino-mediated and a simple gravity-mediated 
model. Moreover, we comment on the case of very heavy lst/2nd generation sfermions. 

^See also the appendix of [3] in this context. 

^For related work in non-minimal supersymmetric models, see also [11, 14-18]. 

^For implications for non-minimal SUSY models with GUT-scale boundary conditions, see [26-30]. 
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Conclusions are contained in Section 5. In three appendices, we give more details about 
the method we use to solve the MSSM renormalization group equations, and comment on 
fine-tuning and on the danger of introducing charge- and color-breaking minima in the 
scalar potential. 

2 Maximal stop mixing and the Higgs mass 

We follow the sign conventions of [32] for soft terms, in which the MSSM superpotential is 
W^Vu UQHu + yd DHdQ + ye EHdL + /x HuHd , (2.1) 
and the SUSY-breaking Lagrangian reads 

>Csoft = - 9 X] + h.C. - (A^yu U QK + Adyd DhdQ + A^ye EhdLj + h.c. 



2 

a=l 



— Bj^ huhd + h.C. + scalar soft mass terms. 

Factorizing the trilinear couplings into Yukawa matrices y and soft coefficient matrices A, 
as in Eq. (2.2), is convenient for our purposes. In fact, for much of what we have to say, the 
only relevant trilinear soft terms are those for the third generation, in particular At = ^-^^33. 
For simplicity we will mostly assume flavor-universal soft terms in the following, in which 
case Au^d.e = ^ol at the GUT scale. 

In the decoupling limit rriA ^ rnz, taking into account the one-loop corrections from 
the top-stop sector, the mass of the lightest MSSM Higgs is 

2 2,tree , a 2,1-loop 



= m^cos^2/3 + log^ + 1 - --^ . (2.3) 



Here as usual tan (3 = Vu/'^di ^ — y^u^ ^d ~ GeV, rrit is the running top mass at the 
scale mt, and M| = mi_^mi^ with m^^ ^ ^^e stop masses. Xt is the stop mixing parameter, 
defined at the scale Ms as 

Xt^ At - iicotp . (2.4) 

The tree-level bound ruh < ruz quickly saturates for tan/3 > 5. To further lift from 
m\ — (91 GeV)^ to around (125 GeV)^, radiative corrections nearly as large as the tree- 
level value are required. In Eq. (2.3) the contribution from the logarithmic term can be 
increased by simply raising M^, but naturalness demands that the soft mass scale should 
be not too far above the electroweak scale. The Xt contribution is easily seen to be 
maximized at \Xt/Ms\ — \/6 = 2.45, although some studies taking two-loop effects into 
account suggest a maximal mixing contribution closer to \Xt/Ms\ ^ 2 (see e.g. [33] for a 
detailed discussion). Values larger than about \/6 will however induce dangerous charge- 
and color-breaking minima in the scalar potential, as detailed in Appendix C. For these 
reasons, in the present paper we will focus on the range 

Ms 



1.5 < 



< 2.5 (2.5) 
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^t^^\yt?At + —^glM, + ... (2.6) 



which we take to define maximal stop mixing. 

If maximal mixing is to be obtained from a GUT-scale model, this places non-trivial 
restrictions on the GUT-scale soft terms. The At parameter generically changes drastically 
during renormalization group (RG) running, because it receives large radiative corrections 
from gluino loops. At one loop, the dominant terms in its RG equation are 

Sizeable gluino masses, which are now favored in the light of direct LHC search bounds, 
will drive At towards large negative values at the electroweak scale (in a phase convention 
where M3 is positive). The soft-breaking masses mg^ and m^^ entering Ms will also receive 
large radiative corrections, because they carry color and because of the large top Yukawa 
coupling. The corresponding one-loop RGEs are 



8^ 



\ybmml^ + ml^+ml^ + \Ai,\^) + ... , (2.7) 

= ^\yt\' i^k + + + 1^*1') - i^almi^ + ... (2.8) 

It is evident that a sizeable M3 will also have a large effect on the scalar soft masses. 

To better quantify the effects of RG running, let us assume universal GUT-scale gaug- 
ino masses Mi/ 2 as predicted by many GUT models, and universal trilinears Aq. The 
most relevant parameters are actually M3, M2 and At, so we are effectively imposing 
M2(Mgut) = Ms{Mqut)'^ We can solve the two-loop RGEs [34] semi- numerically and 
express the electroweak-scale mixing parameter as a function of the GUT-scale soft terms 
(see Appendix A). For Mqut = 2 • 10^^ GeV, Ms = 1 TeV and tan/3 = 20 we obtain 

Xf ^ 9.4 - ^-^ ^0 + 2.2 Al mIi^ - 0.3 AI M1/2 



(2.9) 



(2.10) 



1/2 - u./ ^0 

Terms with coefficients < 0.2 have been suppressed. Similarly, we have 

- 0.4 Al Mlj^ - 0.9 fhl^ + 0.8 ^ . 

The coefficients in these equations vary at most by about 10% in the range 5 < tan /3 < 40. 
The only exception are the coefficients involving /i in Eq. (2.9), which grow as 1/ tan/3 for 
smaller tan/3 as is evident from Eq. (2.4). Explicitly, for the same parameters but with 
tan /3 = 5 we obtain 

^ 9.2 Mi^/2 - ^-^ ^0 Mf/2 + 2.1 M^/^ - 0.3 Al M1/2 

+ 4.3 Mf/2 - 2.5 ^0 M + 0.5 A^ M1/2 /i + 0.8 M^^^ + 0.3 Aq M1/2 , 

(2.11) 



^For simplicity of notation, in the following we will use hatted symbols such as rhn^ , /x, etc. to denote 
GUT-scale boundary values of some of the running MSSM parameters. The exceptions are M1/2 = Mi = 
M2 = M3 and Aq = At = Ab = Ar by definition and, where appropriate, a universal scalar mass mo = 
mo(MGUT). 
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Figure 1. \Xt/Ms\ as a function of with M1/2 the largest GUT-scale soft parameter. 

The upper, violet curve is for tan /3 = 5, while the lower, red curve is for tan j3 = 20. The other soft 
masses are fh'jj^ = m|^^ = = and ttiq^ = m^^ = fn?D^ = 0, in order to maximize Xt in 

Eq. (2.9) while minimizing Ms in Eq. (2.10). Even with this optimal choice of parameters, there is 
no maximal mixing. 



and 



Ml ^ 8.7 + 2.5 M^/^ fhl^ + 1.6 MI/^ + 1.1 MI/^ 
- 0.4 Al Ml,^ - 1.0 m2/2 fhl^ + 0.8 fhl^ . 



3 

(2.12) 



From these expressions one can read off some conditions on the GUT-scale soft terms 
under which maximal mixing results. We will make the extra assumptions that no soft 
term is hierarchically larger than the gaugino mass (such that all terms which we neglected 
in Eqns. (2.9)-(2.12) are indeed subdominant)^ and that all scalar soft masses-squared 
except possibly fn\ju ^^Hd positive.^ It is convenient to classify all possibilities by 
the largest GUT-scale soft parameter among those appearing in Eqns. (2.9)-(2.12). 

• If the largest GUT-scale soft parameter is then maximal mixing is excluded. 
Specifically, for > -M1/2 and (m/f^, m/f^, ffiQ^, fhu^, /i) < M1/2, we find \Xt/Ms\ < 
1.4, as is depicted in Fig. 1. 

• If the sfermion masses are universal, specifically toq^ = fhij^ = tuq, then maximal 
mixing does not allow for mg to be the largest soft parameter. However, maximal 
mixing is possible if either ttiq^ or m^^ is large and the other is small. In that case 
a necessary condition on the spectrum is 

mQ3(f;3) >2Mi/2 (2.13) 



^The relatively large coefficients of the M^/2'^U3 ^i/2^%3 terms in Eqns. (2.10) and (2.12) have 



^Very large scalar soft masses for the first two generations may be an interesting and natural alternative 
scenario for maximal mixing, see e.g. [35-37] and our Section 4.3. 

l/2^L^3 ^1/2^Q3 

led the authors of [38] to propose tachyonic GUT-scale soft masses for the scalars of the third generation. 
This would render Ms small and thus allow for a sizeable \Xt/Ms\ ratio. In this paper however we prefer to 
restrict ourselves to models where sfermion masses are positive at all scales, to avoid possible complications 
from introducing additional vacua to the scalar potential. 
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mus/Mi/2 mus/Mi/2 

Figure 2. Contours of \Xt/Ms\ as a function of mQ^/Mi^2 (vertical axes) and of mu^/Mi^2 
(horizontal axes), for tan/3 = 20 and for Aq = (0, —1, —1.5, —2) x max(mf/,mQ) in the top left, 
top right, bottom left, bottom right panels, respectively. The other soft masses are /i = ttih^ = 
friHd = max(mf/3, mQg). Maximal mixing occurs in the red areas, which shrink to nearly zero as 
\Ao\ is reduced (cf. top left plot). 

and sizeable negative Aq is strongly preferred (even though there remains a tiny 
slice of parameter space where Aq can be zero, if all other parameters are chosen 
optimally). When deviating from the optimal case of sizeable riiH^ and negligible 
^f/sCQs)' required ffiQ^{U3) '• ^1/2 ratio can grow very large. This is illustrated in 
the top row of Fig. 2. 

• 1 can easily be the largest soft parameter if Aq is negative. In the limit that all 
other soft terms are negligible, we find that values of 

Aq ^ -(1-3) max (Ml/2, mQ3,mt/3) (2.14) 

are generic in situations with maximal mixing, particularly if the scalar masses are 
unified or negligible. This is easily understood in the light of the large coefficients 
of the Aq terms in Eqns. (2.9) and (2.10). It is illustrated in Fig. 3 for the case of 
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Ao/Mi/2 



-5 -4-3-2-1 1 

Ao/Mi/2 



Figure 3. Contours of \Xt/Ms\ as a function of A0/M1/2 and of fhH^/Mii2^ with friHd = for 
simplicity and with negligible mfj and ttiq. Left panel: /i = 0, right panel: /i = fhu^' The dashed 
line indicates the border of the electroweak symmetry breaking region according to Eq. (2.16). In 
the red region electroweak symmetry is broken and stop mixing is maximal. 

negligible m^^ and ^gg, and in the bottom row of Fig. 2 for the case of dominant 
sfermion masses friQ^^Us > ^1/2- 

• If fnj^^ is positive, and friHu is the largest GUT-scale soft parameter, then maximal 
mixing seems possible, at first sight, even without significant Aq contributions as can 
be seen from Fig. 3. This is because of the negative-sign M^^fh?^^ contribution in 
Eq. (2.10). However, maximal mixing in this case requires a moderate hierarchy, 

^H„>3Mi/2, (2.15) 

which becomes more pronounced if ffiQ^ and m^/g are non-negligible or if is posi- 
tive, and weaker if Aq is negative. Such a hierarchy is in conflict with electroweak sym- 
metry breaking, as can be understood from the equivalent formula for the electroweak 
symmetry breaking (EWSB) order parameter ruz (here quoted for tan/3 = 20): 

m|^-2 (ImP + ^DIm^ 

= 0.2 - 0.7^0 Ml/2 + 2.9Mi^/2 " 2-1 \U? - ^-^rhj^^ + 0.7 + 0.8 

(2.16) 

In RGE language, the gaugino masses (specifically the gluino mass) are typically 
responsible for driving Toj^^ negative at the electroweak scale, thus triggering elec- 
troweak symmetry breaking. If the GUT-scale value for Toj^^ is too large, and M1/2 
is too small, electroweak symmetry is not broken, because the RHS of Eq. (2.16) 
remains negative. Sizeable friQ^ or m^/g, or sizeable negative Aq can remedy this, but 
only the latter is favorable for maximal mixing. 

Figs. 2 and 3 are based on tan/3 = 20, cf. Eqns. (2.9) and (2.10). The situation remains 
the same qualitatively also for smaller tan/3. For tan/3 ^ 5, relatively large /i can however 
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Figure 4. Contours of \Xt/Ms\ as a function of /i/Mi/2 and tan/3, for friH^ = f^Hd = ^1/2 and 
mQg = mf/3 = 0. Left panel: A0/M1/2 = — 1, right panel: A0/M1/2 = —3/2. Clearly, at small 
tan /3, large /i can lead to maximal mixing in regions which are otherwise at the borderline. 

slightly widen the allowed regions for maximal mixing, since Xt — At — ji/ tan/3 has a 
stronger //-dependence if tan/3 is small. The effect of /i on maximal mixing at smaller 
tan/3 is also illustrated in Fig. 4. We do not consider values of tan/3 < 5, since they no 
longer maximize the tree-level Higgs mass Eq. (2.3). 

To summarize, generically, maximal mixing in a GUT-scale model with unified gaugino 
masses at tan /3 > 5 requires 

large negative ^ —(1-3) max (Mi/25 ^Qs? ^t/s)- 

This is a non-trivial requirement on any GUT-scale model. In addition, it is beneficial but 
not strictly necessary for maximal mixing if there are 

positive up-type Higgs soft masses fh?^^, 

and 

small third- generation soft masses fnfj^, ^Qs' 

as compared to the gaugino mass. 

The resulting Xt/Ms at the electroweak scale will necessarily be negative. All this we 
have deduced from the semi- numerical solution of the RGEs underlying Eqns. (2.9), (2.10), 
(2.11), and (2.12); see Appendix A for more details. The picture is confirmed by parameter 
space scans performed with SoftSusy3.2.4, whose results are shown in Section 4. 
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3 Models 



We assume F-term SUSY breaking in some hidden sector, mediated to the visible sector 
by messenger states which can be supersymmetricahy integrated out at some high scale 
M. If X = FO'^ is the Goldstino background field which parametrizes SUSY breaking, the 
lowest-order operator inducing a gaugino mass takes the form 

Cm,,, = ci/2 j d^e^ tr W^^W^, + h.c. (3.1) 

As before we assume gaugino mass unification, i.e. a GUT-preserving F-term VEV. 

If the hidden sector couples to the Higgs, we can have the Giudice-Masiero operator 

C^^c^ J d^e ^HuHd + h.c. (3.2) 
which induces a ji term, and the operators 

= j d^e^ (C^„ + |^^|2) (3 3) 

which contribute to the Higgs soft masses. The operators 

Ca^ j d^e^ (c^„ \Hu\' + CA, m') + h.c. (3.4) 

are often neglected because they can be absorbed by a holomorphic field redefinition 
Hu^d Hu^d + / ^) ' order to cleanly separate hidden and visible sectors, 

we instead keep these operators explicit, noticing that they also induce soft Higgs masses, 
and (together with the Yukawa couplings) flavor-universal trilinear A-terms, and (together 
with Cjj,) a Bf^ term. If ca^ and ca^ were set to zero by the above field redefinition, these 
soft terms would instead arise from the induced superpotential terms and the change in 
the ch^^ coefficients in Eq. (3.3). Finally, there are the operators 



/X^ X 
^^^J^HuHd+ h.c. (3.5) 



which also induce a term. 

If the hidden sector couples to matter fields, the equivalents of Eq. (3.3) and (3.4) 
with Higgs fields replaced by matter fields can be present. They will induce soft masses and 
flavor non-universal A-terms. The latter can also arise from superpotential operators of 
the form 

/ d^e— {QHuU + QHdD + LHdE) + h.c. (3.6) 

We will assume that the ^/S^ problem is solved; in particular there are no unacceptably 
large contributions to /i and as would arise from the renormalizable terms 



/ 



(f9 iflHuHd + XHuHd) + h.c. (3.7) 
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From the discussion in Section 2 it is clear that, in order to reahze maximal mixing, 
models with vanishing or suppressed direct couplings to matter fields are preferred. The 
operators of Eqns. (3.1)-(3.5) are then all that is needed to fix the high-scale soft masses. 
Furthermore, the model should allow for Ci/2 ^Au be chosen such that the resulting 
At : Ml/2 ratio is around —(1-3). SUSY breaking scenarios with suppressed couplings of 
the hidden sector to matter fields are not generic, because the scalar masses cannot easily 
be forbidden by symmetry (by contrast, suppressed gaugino masses are easily obtained 
but phenomenologically undesirable). We will list a number of well-known examples, and 
comment on the implications of our study. 

3.1 Gaugino (-Higgs) mediation 

Gaugino-mediated supersymmetry breaking in its minimal form is defined by the gaugino 
masses being the only non-vanishing terms at the mediation scale. In other words, only the 
operators of Eq. (3.1) are present at the scale M; the only MSSM fields to couple directly 
to the hidden sector are the gauginos. This was originally motivated by 5D models [39-41] 
in which the hidden sector and the chiral (Higgs and matter) superfields of the MSSM 
were separated in an extra dimension, with only the MSSM gauge fields coupling to both. 
Deconstructed models [42, 43], models with Seiberg duality [44], and models with strongly 
coupled near-conformal hidden sectors (see below) may also give rise to gaugino mediation. 

Independently of maximal mixing, the minimal scenario is not viable phenomenolog- 
ically because it is missing a jn term. Realistic extensions therefore require that also the 
Higgs fields should be directly coupled to the hidden sector. In the 5D picture, the Higgs 
and gauge fields would be bulk fields, while the MSSM matter fields would be localized on 
one brane and the hidden sector on the other. Generically, all operators in Eqns. (3.1)- (3.5) 
are then present at the scale M, and maximal mixing can be naturally accommodated. In 
fact, all three conditions listed at the end of Section 2 are satisfied if ca^ — — 2ci/2- 

It should be mentioned that the mediation scale M can be parametrically lower than 
Mgut in gaugino-mediated models. For significantly lower M, running effects become less 
important, and the condition for maximal mixing eventually approaches Eq. (2.5) with 
the boundary values of At and fi at the scale M substituted, and with mq^ and m^/g well 
approximated at leading-log order. 

3.2 Scalar sequestering 

Scalar sequestering [45, 46] is a more restrictive version of gaugino-Higgs mediation. The 
hidden sector is assumed to be close to a strongly-coupled conformal fixed point over a 
large range of energies, with X a composite operator satisfying A^t^ — 2Ax > 0. This 
implies that [46, 47]^, at the scale M where the theory exits the strong coupling regime, 

Cl/2 ^ CA^,d ^ > (cH^,d - l^A^J^ - |C^P) - (CB^ - C^{ca^ + C^J) (3.8) 

The operators responsible for matter soft masses also end up being suppressed. Eventually 
the dominant soft terms at the scale M are M1/2, Aq, fh?^ ^ and /i, which are all of the 

^In [48] it was shown that, strictly speaking, these relations do not hold generically but depend on 
additional assumptions about the hidden sector. 
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same order, while the combination 

;2 , ^2 



m 



k.. +W -B, (3.9) 

and the matter soft terms are suppressed. 

Our analysis is applicable for M close to Mqut- The phenomenological consequences 
of the condition Eq. (3.9) were discussed in detail in [49]. It was shown that for > 0, 
almost always has the same sign as At (and the opposite one if < 0). Moreover, small or 
vanishing requires also small At, in apparent conflict with maximal mixing, cf. Fig. 3 
in [49]. 

3.3 Radion mediation / Scherk-Schwarz SUSY breaking 

In 5D models, supersymmetry can be broken by giving an F-term to the radion multiplet. 
This multiplet hosts the 5D gravitational degree of freedom which corresponds to the 
compactification radius. Radion mediation is equivalent to breaking SUSY by the Scherk- 
Schwarz mechanism, and generalizes to modulus-mediated SUSY breaking in superstring 
models. Of particular interest are models whose compactification scale is close to Mqut, 
since in that case boundary conditions can be used to break the grand-unified symmetry 
down to the MSSM. Minimal models of this kind give rise to specific GUT-scale soft term 
patterns, so it is natural to ask if they can also accommodate maximal mixing. 

The role of the operator X in Eqns. (3.1) to (3.5) is then played by X = TM/{2R), 
where M is identified with the 5D Planck mass, and T is the radion multiplet with (T) = 
i? + F^9^. For the simplest model with a flat 8^/1^2 extra dimension, the gaugino masses 
are (see e.g. [50]) 

Mi/2 = |^, (3.10) 

and the A-terms depend on the localizations of the matter fields and on the origin of the 
Higgs field. They are given by the sum of the contributions from the Higgs and the matter 
fields involved in the respective trilinear coupling. Roughly speaking, in gauge-Higgs unified 
models where the Higgs originates from the 5D gauge multiplet, the Higgs contributes 

aA = (3.11) 

but bulk matter fields Q3 and U-^ originating from 5D hypermultiplets each contribute 

= +^ (3.12) 

leading to At — i^i the most naive model with unlocalized and ^73. Localizing 

and U-^ towards one of the branes brane allows to reduce their contribution, but in 
the potentially interesting limit in which they are completely brane-localized (which would 
leave us with At = —F^ = —^1/2 from the Higgs contributions), the Yukawa couplings 
vanish. In models without gauge-Higgs unification, with the Higgs coming from a bulk 
hypermultiplet, it gives a wrong-sign contribution 

Alt = +^. (3.13) 
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In short, maximal mixing does not occur in minimal radion-mediated models. While in 
realistic models the exact sfermion soft terms depend on the modelling of the matter sector, 
and more elaborate examples may allow for maximal mixing in principle, we were unable 
to find an example. For instance, for the model of [51] we find that is bounded 

such that it can never become large and negative, ^ —0.72. Furthermore, the 

third-generation soft masses are typically comparable and of the order of the gaugino mass, 
so large ratios which might still allow for maximal mixing (see the upper row of Fig. 2) do 
not appear. In conclusion, maximal mixing seems not to be a generic feature of radion- 
mediated SUSY breaking in 5D. 

3.4 Gauge mediation 

Gauge-mediated supersymmmetry breaking, by definition, encompasses models whose hid- 
den sector decouples from the MSSM as the MSSM gauge couplings are switched off [52]. 
Similar as in gaugino- mediated models (with which there is indeed some overlap), fi is 
missing in pure gauge mediation, and needs to be generated by additional Higgs-hidden 
sector couplings. This generically induces a too large S^, a problem which needs to be 
solved in realistic models, or overcome with very special soft mass patterns [53, 54]. 

Pure gauge-mediated models also predict vanishing A-terms at leading order at the 
mediation scale. A priori it is not clear that this prediction is maintained when the model is 
extended such as to solve the /i/S^ problem, but it was shown in [55] that generically this 
is indeed the case (although Higgs-messenger couplings in the superpotential can be used 
to obtain large negative A-terms while leaving the /i/S^ problem unsolved [56]). Therefore 
models of gauge-mediated supersymmetry breaking do not lead to maximal mixing. 

The mediation scale M is often taken to be far below Mqut in gauge-mediated models 
(for exceptions, see e.g. [57, 58]). However, this generically does not improve the prospects 
for achieving maximal mixing (see also [9]). The implications of a 124-126 GeV CP-even 
Higgs boson for minimal gauge mediation have been investigated in detail in [59] , one of the 
conclusions being that "the majority of the sparticle masses are in the several to multi-TeV 
range" . 
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4 Numerical analysis 



Let us now illustrate the impact of maximal mixing by means of parameter space scans 
of two models. The first is gaugino mediation, and the second is a more generic gravity- 
mediated model with non-universal Higgs masses (NUHM) and universal sfermion masses 
mo. In the NUHM case, we will distinguish between mo < ^ind mo > Mi gaugino 

mediation can be regarded as a special NUHM scenario with mo = 0. In addition, we 
will also briefly survey a NUHM scenario where maximal mixing is generated from large 
sfermion masses for the first two generations. 

We use SoftSusy3.2.4 [60] for the spectrum calculation, micrOMEGAs2.4 [61] for 
low-energy observables, and Hdecay4.4 [62] for computing Higgs decays. Besides radiative 
EWSB and the absence of tachyons, we require the following phenomenological constraints 
to be satisfied: 

• Mass limits: m~± > 103 GeV, ruf^ > 92 GeV, mg^^ > 100 GeV, m^^ > 100 GeV, 
rUg > 500 GeV, and rrih > 115 GeV to avoid the most stringent constraints from 
collider searches, as well as 

• Flavor constraints: 2.87 x 10"^ < BR{B < 4.23 x 10"^ [63] and BR{Bs 
IJi^ lJi~) < 5.4 X 10~^ (i.e. the LHCb limit [64] augmented by a 20% theoretical uncer- 
tainty) to be compatible with recent S-physics results. 

We furthermore compute the dark matter relic density Q.h^ and the SUSY contribution 
Aa^ to the muon anomalous magnetic moment, but do not impose any restrictions on them 
a priori. 

We will also discuss the implications of maximal mixing for a possible Higgs signal 
near a mass of 125 GeV. To this end, we approximate the signal strength for a given final 
state X, relative to the Standard Model expectation for the same Higgs mass, as 

o/J^^ ^ [<y{99^h) BR(/i ^ X)]mssm _ [r(/i ^ gg) BR(/i ^ X)]mssm u n 

^ ^ [a{gg^h)BR{h^X%^ ^ [V{h ^ gg)BR{h ^ ' ^'^ 

This is justified because differences in (j{gg h) versus V{h gg) should largely cancel 
out when taking the taking the MSSM/SM ratio. It is important to note that SUSY 
contributions can lead to modifications in both Higgs production and Higgs decays as 
compared to the SM. The effective ggh coupling is dominated by the top-quark loop, 
while the effective /i77 coupling is dominated by the contribution from W bosons with a 
subdominant contribution of the opposite sign from top quarks. Both couplings can receive 
a large contribution from third-generation sfermions, in particular from stops [66]. In case 
of no stop mixing, the light stop loop interferes constructively with the top loop, while 
the interference is destructive in the case of large mixing. In the first case hgg rate will 
increase, while in the latter case it will decrease. The contrary is true for the h^^ couplings. 
Moreover, weakly interacting particles such as charginos and sleptons will contribute only 
to the /i77 coupling. In particular, light staus that are strongly mixed can enhance the 
/i ^ 77 rate [10]. Another important effect is the enhancement or suppression oi h ^ bb, 
which can significantly change the overall branching ratios. In fact an enhancement of the 
77 signal, i.e. ^^(77) > 1, is often due to a suppression of BR(/i bb). 
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Figure 5. Projected parameter space in the gaugino mediation model. The amount of stop mixing, 
Xt/Ms^ is shown by a color code. Brighter (darker) shaded points of each color have rrih > 123 
{rrih < 123) GeV. Evidently, maximal mixing requires A0/M1/2 — —1 to —3. 

4.1 Gaugino mediation 

In the gaugino- mediated model the non- vanishing soft masses at the GUT scale are 
Aq, /i, 5^, Tnj^^ and irij^^- Given m^, the last four parameters can be traded against 
tan/3, /X, and the pseudoscalar mass tua- Our input parameters are therefore M1/2 and Aq 
at Mgut, together with tan/3, fi and ttia at the electroweak scale. We choose fi > and 
perform flat random scans letting M1/2, M vary up to 2 TeV; tan /3 is allowed to vary 

between 1 and 60. Regarding Aq, we scan over two different intervals, ^ 3 and 

Aq G [—4, 0] TeV. Of 256k valid points from these scans, 158k remain after the basic mass 
limits and the flavor constraints listed above. Of these, 10k points have rrih = 123-127 GeV. 

Figure 5 shows a projection of the scanned parameter space in the Aq versus M1/2 
plane, with the amount of stop mixing visualized by a color code. For each color, marking 
a certain interval of stop mixing, points with rrih > 123 (rrih < 123) GeV are shown in 
brighter (darker) shades. Evidently, maximal mixing (yellow, red and blue points) requires 
Aq/Mi/2 — —1 to —3. As can also be seen, a SM-like h with mass ruh > 123 GeV requires 
large mixing with a large negative Aq, or very large M1/2. Even with \Xt/Ms\ ^ 2, a Higgs 
near 125 GeV requires M1/2 > 750 GeV; for \Xt/Ms\ 1, M1/2 is pushed up to 2 TeV. 

The correlation of ruh with the stop mass scale and mixing is further illustrated in 
Fig. 6. We see that light stops with Ms ^ 1.5 TeV indeed require maximal mixing, 
—2.5 < Xt/Ms ^ —1.5, for having a Higgs mass compatible with the ATLAS and CMS 
excesses. We also see that while in principle it is easy to have Ms < 1 TeV as preferred 
by naturalness, once we impose rrih > 123 GeV the stop mass scale is pushed to a TeV 
and above; if we demand Ms ^ 1 TeV in addition to rrih = 123-127 GeV, this requires 
-2.5 < Xt/Ms < -2, see the RHS plot in Fig. 6. 

The ruh dependence on various quantities (tan/3, Mi/2^ m^^ and Aq/Mi/2) is shown 
explicitly in Fig. 7. Maximal mixing clearly prefers not too large tan/3 < 25, while 
ruh > 123 GeV needs M1/2 ^ 750 GeV. The highest h mass is achieved for Aq/Mi/2 ^ —1-5, 
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Figure 6. Higgs mass rrih versus stop mixing Xf/Ms (left) and Xt/Ms versus Ms (right) in the 
gaugino mediation model. 
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Figure 7. Lightest Higgs mass dependence on various quantities, in the gaugino mediation model. 
Top row: Large rrih favors not too large tan/3, especially when demanding a lighter spectrum, and 
^1/2 ^ 750 GeV (i.e. rrig ^1.7 TeV). Bottom row: ii masses below 1 TeV can still be reconciled 
with large m^, but only at maximal mixing. If the gluino is light enough to be seen at the LHC, 
maximal mixing is also strongly favored. 
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Figure 8. Scatter plots of points with rrih = 123-127 GeV in the gaugino mediation model, on 
the left in the ii versus g mass plane, on the right in the q versus g mass plane. Since mo = 0, 
rriq < rrig throughout the parameter space. 

but stops are heavy in this case, above 2 TeV. Although large rrih still allows for ii masses 
below 1 TeV in case of maximal mixing, overall a Higgs in the 123-127 GeV mass range 
points towards a heavy spectrum. The correlations between gluino, stop and l^^/2^^ 
generation squark masses are shown in Fig. 8. Remarkably, the latest LHC bounds of 
rug^q ^1.4 TeV [67, 68] for CMSSM-like scenarios with rrig 2^ rriq are automatically evaded. 
In fact, requiring ruh > 123 GeV, we find m^^ > 715 GeV, ruq > 1.5 TeV and rUg > 1.7 TeV 
in gaugino mediation. Moreover, for tan/3 < 50 we find a maximal h mass of around 
125 GeV (which is however subject to a 1-2 GeV theoretical uncertainty). 

Owing to the vanishing sfermion soft masses, mo = 0, over most of the gaugino medi- 
ation parameter space, the LSP is the lighter stau fi.^ (92% of the points satisfying mass 
limits, flavor constraints and rrih = 123-127 GeV have a fi LSP, while 7% have a neutralino 
LSP.) Hence one might expect that the gg ^ h ^ rate be enhanced by the light fi 
contribution [10]. While the /i ^ 77 partial width is indeed enhanced for light staus^ (and 
light stops with large mixing), this is mostly compensated by the suppression of the h ^ gg 
partial width due to the i loop contribution, and by the larger total h width (mostly due 
to larger T(h bb)). As a result, the gg ^ h ^ signal strength is typically 80-90% 
of that in the SM. A similar suppression arises for the ZZ final state. This is illustrated 
in Figs. 9 and 10. The few points with a signal strength i? > 1 at rrih > 123 GeV feature 
very large tan /3 = 52-60, and a reduced h ^ bb rate. 

Finally, in Fig. 11, we illustrate the implications of maximal mixing and a heavy CP- 
even Higgs for the GUT-scale Higgs soft masses m/f^, m/f^, and for the weak-scale Higgs 
mass parameters /x and rriA- (We use the convention mn^ ^ = sign(m|^^^) x ^\rri\j 

As can be seen, positive up-type Higgs soft masses fn\ju preferred, in particular in case 
of small fji (as preferred by fine-tuning). This confirms our expectations from section 2. 

^More precisely, to obtain a consistent cosmological picture the stau should in that case be the next-to- 
LSP, and the true LSP a gravitino or axino [65] . 

^Note that rrih > 123 GeV leads to rrif-^ > 99 GeV in our dataset. 
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Figure 9. Partial widths of h ^ 77, h ^ bb and h gg, relative to SM expectations. 




Figure 10. In the gaugino mediation model with maximal mixing, the gg ^ h ^ rate is slightly 
below the Standard Model one (left), mostly as a result of the reduced gluon fusion production cross 
section. On the right, we show the ZZ signal strength. The few points with a signal strength R > 1 
at rrih > 123 GeV have very large tan/3 . 

Moreover, the pseudoscalar Higgs mass ruA prefers to be large in the maximal mixing case, 
well above current limits. 
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Figure 11. Correlations between GUT-scale Higgs soft masses friH^ and m//^, and weak-scale 
Higgs mass parameters /i and m^. 
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Figure 12. Projected parameter space in the NUHM model, analogous to Fig. 5. We distinguish 
Ml/2 > ^0 (left) and mo > M1/2 (right). The amount of stop mixing, X^/M^, is indicated by a 
color code, with brighter (darker) points of each color having > (<) 123 GeV. It is apparent 
that negative Aq is necessary for maximal mixing, with the ratio Ao/max(Mi/2 7 ^0) ^ —1 to —3. 

4.2 NUHM model 

Let us now turn to the gravity-mediated model with non-universal Higgs soft masses 
(NUHM). The parameters are the same as in the gaugino- mediation case, but for non- 
vanishing soft masses for squarks and sleptons at Mqut- For simplicity, we take the latter 
two to be universal at the GUT scale. (As discussed in section 2, the scalar masses which 
actually affect maximal mixing are mainly tuq^ and mu^.) We scan over ^o, tan/3, 

/i and rriA as in the previous subsection, allowing however l^ol up to 6 TeV. In addition, 
we let mo vary from to 5 TeV. Over most of the parameter space, the non-vanishing mo 
makes the leptons heavier than the Xi- We thus keep only points with a neutralino LSP, 
without however restricting fi/i^. We will discuss the Xi relic density at the end of this 
subsection. 

Figure 12 shows a projection of the scanned NUHM parameter space, analogous to 
Fig. 5 in the gaugino mediation model. As discussed in section 2, it makes a difference 
whether M1/2 or mo is the largest soft mass. In Fig. 12 and following figures, we hence 
distinguish between the two cases M1/2 > ^0 and mo > Mi/2- As expected, maximal 
mixing requires a ratio between Aq and max(Mi/25^o) of about —1 to —3. Moreover, 
a heavy (rrih > 123 GeV) MSSM Higgs requires large mixing with a large negative Aq 
— the tip of the scatter plot being again around Aq ^ —2 TeV — or an overall heavy 
spectrum. A difference to the case with vanishing mo is that now much larger negative 
values of Aq still give a valid spectrum. Another interesting difference is that with increasing 
mo, larger values of \Xt/Ms\ become consistent with a heavy h. Indeed, for large mo, 
the lowest Ms giving ruh = 123-127 GeV is found for Xt/Ms ^ -3 to -3.5. On the 
other hand, as previously mentioned such large values of \Xt/Ms\ give rise to dangerous 
charge- and color-breaking minima in the scalar potential, so although SOFTSUSY gives 
a valid spectrum the corresponding points should not be trusted to be phenomenologically 
viable. The correlations between m/^, Ms and amount of stop mixing are shown in Fig. 13. 
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Figure 13. Higgs mass rrih versus stop mixing Xf/Ms (top row) and Xf/Ms versus Ms (bottom 
row) in the NUHM model. 

Furthermore, Fig. 14 shows the dependence of rrih on tan^ and m^^, with Xt/Ms again 
indicated by a color code. We see that at large mo, maximal mixing is possible also for 
large tan /3. As a side remark we note that the maximal h mass in these scans is about 
128 GeV, consistent with the current 95% CL limit. 

Consequences for LHC SUSY searches are illustrated in Fig. 15. For M1/2 > ^o, we 
still find ruq 2^ ruq ^1.5 TeV. For mo > however, gluinos can be as light as 500- 

600 GeV (with stops being light and maximally mixed, and rrih iii the desired range). First 
and second generation squarks need to be heavy in this case, around 2-3 TeV, as can be 
seen in the bottom-right panel in Fig. 15. 

Expectations for the Higgs signal in the 77 and ZZ channels are shown in Fig. 16. 
We observe that for a Higgs mass in the desired range, i?(77) and R{ZZ) < 0.9, rather 
independently of the stop mass (a decoupling effect of heavy stops can however be seen 
in the lower boundary of i? at a given Xt/Ms)- Higgs signal strengths close to or above 
1 occur for mo ^ ^1/2 5 they require heavy stops with small mixing, combined with large 
tan/3 and large m^, such that the h ^ bb rate is suppressed. 

Finally, although it does not directly have to do with maximal mixing, let us consider 
the question of neutralino dark matter. The relic density of the neutralino LSP is plotted 
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Figure 14. Dependence of rrih on tan (3 (top row) and on m^^ (bottom row) in the NUHM model. 
As mo increases, larger tan (3 values are consistent with maximal mixing. Moreover, ii masses 
below 1 TeV can easily be reconciled with large rrih ~ 125 GeV if mo and the stop mixing are large 
enough. 

versus the neutralino mass in Fig. 17. Interestingly, for M1/2 > a large fraction (45%) 
of the points have Qh^ < 0.135, and overall the relic density does not exceed 20. The 
points with very low fi/i^ typically feature a higgsino-like LSP, which makes the scenario 
difficult to detect at the LHC [69, 70]. In the remaining cases, when fi is large, Qh^ is 
low because of co-annihilations. For mo > M1/2, the situation is quite different, and we 
find the "usual" MSSM picture with nh'^ ranging from 10"^ to 10^. Roughly 2% of the 
points have 0.09 < fth'^ < 0.135. An example of a "perfect" point the sense of light stops, 
maximal mixing, rrih — 125 GeV and O^o/i^ = 0.1 is given in Table 1. 
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Figure 15. Scatter plots of points with ruh = 123-127 GeV in the NUHM model, in the top row 

ii versus g masses, in the bottom row q versus g masses. 
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Figure 16. In the NUHM model with maximal mixing and rrih = 123-127 GeV, the gg ^ h ^ 
and ZZ rates are only about 60-90% of those in the Standard Model one, mostly as a result of the 
reduced gluon fusion production cross section. An enhancement {R > 1) only occurs for small stop 
mixing and when the h ^ bb decay mode is suppressed at large tan/3 and large ttia- • 
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Figure 17. Dark matter relic density versus neutralino LSP mass in the NUHM model, on the left 
for Ml/2 > 7 on the right for mo > Mi/2- All points have ruh = 123-127 GeV. 
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Figure 18. In the case of very heavy first and second generation sfermions (here moi = 10 TeV) 
and light stops, sbottoms and staus (mos < M1/2 < 2 TeV), Ms is driven to smaller values and 
maximal mixing occurs even for l^ol < Mi/2- We find rrih > 123 (124) GeV with m^^ < 1 TeV for 
ratios as low as Ao/Mi^2 ^ —0.2 (—0.6). 

4.3 Split generations, inverted sfermion-mass hierarchy 

So far we have considered only the case that there is no large hierarchy between the three 
generations of sfermions. An interesting alternative, motivated also by the absence of any 
signal of new physics in the flavor sector, is the case of an inverted sfermion-mass hierarchy, 
with squarks and sleptons of the first two generations being very heavy (in the multi-TeV 
range) while the third generation and the gauginos are light, of the order of 1 TeV. Together 
with the requirement of small /x, this is often referred to as "effective SUSY" or "natural 
SUSY" in the literature. 

To illustrate this case, we perform a scan over the NUHM parameter space as before, 
but setting the soft masses of the first two generations moi = 10 TeV. For the third 
generation, we assume a universal soft mass mo3, which we let vary between and Mi/2- 
As also shown in [35], in this setup a Higgs near 125 GeV with light stops is possible even 
for small \Ao\. This is because, during RG running, Ms is driven down by the first two 
generation squarks being very heavy [71]. Therefore maximal mixing now occurs also at 
smaller Aq, see Fig. 18. Scenarios with maximal mixing, m^^ below 1 TeV, and a Higgs 
near 125 GeV can now be found for Aq/Mi/2 ratios of around —0.2 to —2. 
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5 Discussion and conclusions 



If the MSSM is to accommodate a Higgs boson around 125 GeV, maximal stop mixing is the 
only way to avoid multi-TeV stop masses. Since the Higgs sector and the top/stop sector are 
coupled to each other by the large top Yukawa coupling, multi-TeV stop masses generically 
imply multi-TeV mass parameters in the Higgs potential. A relatively small electroweak 
scale can then only arise through delicate cancellations, which requires considerable fine- 
tuning. In the context of phenomenological models which prescribe the MSSM parameters 
at the TeV scale, it has therefore been argued that, in the most natural remaining regions 
of the MSSM parameter space, the stop masses should be as small as possible (see, for 
instance, [72, 73]). To be compatible with a 125 GeV Higgs, the stops will then have to be 
maximally mixed. 

However, as we have argued, one of the most appealing aspects of the MSSM is that 
it can be valid up to very high energies. It is not clear a priori if a given set of TeV- 
scale MSSM parameters can result from some healthy UV completion at a high scale, or 
if instead it is a point in the "swampland". For example, and of immediate relevance to 
maximal mixing, it is not possible to obtain Xt/Ms ^ +2 from a GUT-scale model (barring 
hierarchically large GUT-scale trilinear terms), because of the negative gluino contribution 
to At during its RG evolution. Indeed, we have shown in some detail that in models where 
all GUT-scale soft parameters are of the order of the gaugino mass or smaller, the trilinear 
coupling must be large and negative at Mqut foi" maximal mixing to result. In that case 
one has Xt/Ms ^ —2 at the electroweak scale. 

In models of high-scale SUSY breaking mediation, the fine-tuning required to obtain 
a small electroweak scale from a large soft mass scale also involves the gluino mass and 
the top trilinear. The electroweak scale is very sensitive to the UV-scale M3, since gluino 
loops strongly affect the stop masses and thus indirectly the Higgs mass parameters and 
the electroweak scale. This is evident e.g. from the large coefficient of M3 in Eq. (A. 2). 
The sensitivity of the electroweak scale with respect to At is less pronounced. Overall, the 
least fine-tuned remaining parameter regions of the MSSM are characterized by low M3 
(or low Ml 12 if the gaugino masses are universal), correspondingly low Ms^ and maximal 
mixing; see Appendix B. Even in these regions the fine-tuning is at the permille level or 
worse [5, 74]. 

Maximal mixing does not single out any particular scenario for SUSY breaking me- 
diation; it follows from a parameter choice in models where the UV-scale trilinear soft 
terms are free parameters. It would be highly interesting to identify models which actually 
predict large and negative A-terms. Among the classes of models we studied, only gauge 
mediation and to some extent radion mediation predict the trilinear terms, and in these 
cases the prediction disfavors maximal mixing. 

In models where maximal mixing is allowed, it becomes an interesting question what 
its phenomenological consequences are. We have studied three examples, with particular 
attention to the parameter regions which give a Higgs boson around 125 GeV. In each case 
all experimental constraints can be satisfied for suitable parameter choices. Squarks and 
gluinos typically turn out to be heavy, beyond the reach of the 2011 LHC run at ^/s — 7 
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TeV. They may however be within reach of the ^/s = 8 TeV run in favorable cases. The 
gg ^ h ^ signal strength consistently turns out to be some 10-40% below that of the 
Standard Model. 

Table 1 lists the properties of four representative points, two for gaugino mediation 
and two for the NUHM case, with stops below 1 TeV and maximal-mixing. GM-1 is the 
point with lowest M^-, while GM-2 is the point with lowest Ms and < 500 GeV from 
the gaugino mediation scan. Both these points have a stau as the lightest sparticle with a 
relic abundance of the order of 10~^.^^ NUHM-1 is a low-stop-mass point from the NUHM 
scan with M1/2 > rriQ. It features a neutralino LSP with large higgsino component and 
a relic density which is too low, so that the Xi would provide only about 20-30% of the 
dark matter. NUHM-2 has a large rriQ of order 2 TeV and Aq ^ — 2mo. It exactly matches 
the desired Higgs mass (125 GeV) and has a bino-like neutralino LSP with a relic density 
of nh^ = 0.1. The SLHA files of these maximal mixing benchmark points are included as 
ancillary files in this preprint. 
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GM-1 


GM-2 


NUHM-1 


NUHM-2 


mo 








759 


2295 


Ml/2 


772 


974 


833 


736 




-2254 


-2099 


-2336 


-4140 


tan/? 


15.3 


10.1 


14.7 


19.4 


/i 


1343 


273 


329 


1049 


rriA 


1981 


1640 


330 


678 


rrih 


123 


123.4 


123.6 


125.2 


Al 


326 


263 


302 


319 


m~o 


620 


282 


338 


611 




620 


274 


327 


611 




221 


116 


604 


2064 




425 


255 


685 


2210 




469 


675 


992 


2387 


niq 


1532 


1894 


1808 


2681 


m~g 


1712 


2120 


1863 


1749 




730 


821 


652 


862 




1306 


1549 


1396 


1894 


Xt/Ms 


-2.32 


-2.13 


-2.36 


-2.25 




1498 


1521 


-273 


566 




585 


1934 


1991 


3030 


BRb^,^ [10-4] 


2.94 


2.88 


3.16 


3.18 


BRb^^^ [10-9] 


3.14 


3.07 


4.60 


4.26 


Aa^ [10-10] 


6.87 


2.64 


3.47 


1.07 


!^/i2('LSP') 


8.5 X 10-3 


6.5 X 10-3 


0.026 


0.101 


<7si(XiP) [Pb] 






2.9 X 10-^ 


3.7 X 10-10 


<7SD (XiP) [pb] 






1.5 X 10-4 


1.4 X lO-'^ 



Table 1. Sample points with light stops, maximal mixing, and a Higgs near 125 GeV in the 
gaugino mediation (GM-1, GM-2) and NUHM (NUHM-1, NUHM-2) models. ^^/i^CLSP') is the 
relic abundance of the fi for the GM points, and of the Xi for the NUHM points, asi and ctsd 
are the spin-independent and spin-dependent scattering cross sections off protons; for NUHM-1 
these should be rescaled by a factor ^ = l]/i^/0.1123 for comparison with the experimental limits 
(Ccrsi = 6.7 X 10-^ and ^ctsd = 3.5 x 10"^ for NUHM-1). 
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Appendix 

A Semi-numerical solutions of the MSSM renormalization group equa- 
tions 

The one-loop RGES of the MSSM can be integrated analyticahy when keeping only the top 
Yukawa coupling nonzero [75].^^ The procedure can be summarized as follows: As a first 
step, one writes down the closed- form solutions for the gauge coupling and gaugino mass 
RGEs, which are of course easily found and well known. The top Yukawa RGE is a Bernoulli 
equation whose solution, given the solutions for the gauge couplings, can be expressed as 
a simple integral. Then the At RGE becomes a linear ODE with known inhomogeneous 
term and known coefficient functions, and is easily solved in terms of an integrating factor. 
Finally, having solved all of these RGEs, the coupled RGEs for third-generation scalar 
masses and for Toj^^ can be integrated in a suitable basis. 

In this paper we are using a somewhat refined approach to improve precision. We 
keep all third-generation Yukawa couplings, and we are using two-loop RGEs for the gauge 
couplings, Yukawa couplings, and gaugino masses. Boundary values for the gauge and 
Yukawa couplings are matched to SOFTSUSY GUT-scale values, in order to properly take 
threshold corrections into account. With the more complicated coupled two-loop system to 
solve, the solutions can no longer be expressed by simple integrals, but this is unnecessary 
to extract the information we need. We will now briefly describe our method. 

As a first step, we fix some value of tan/3 and some soft mass scale Ms- Using 
appropriate boundary values for the gauge and Yukawa couplings, we solve their RGEs 
numerically between Mqut cind Ms- What we are eventually interested are however the 
SUSY breaking mass parameters. Quite generally, from the structure of the RGEs and 
from dimensional analysis, it follows that their values at Ms take the form 

Ma{Ms) = 5^ «a6 M5 + ^ a'^^A^ , 

b X 

A^ciMs) = J2 ^^oMa + J2 f^-y^y ' (A.l) 

a y 

ab xy xa x 

Here hatted quantities denote boundary values at Mqut as in the main text, and the /3, 7 
coefficients are functions of tan/3 and of Ms- They are obtained by numerically solving 
the RGEs with special boundary conditions. For instance, setting all GUT-scale masses 
except Ml to zero allows to read off the aai, /3xi, and 7(^ii coefficients from the numerical 
solution, and similarly for the others. The same method is also applied to the fi and Bj^ 
RGEs. 

The coefficients we find roughly agree with the one-loop running, tree-level matching 
value often found in the literature, with the exception of M3 whose influence on the low- 
scale soft terms we find to be somewhat reduced. Compare for instance our result for mz 

^^See e.g. [76, 77] for reviews and applications. 
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at tan /3 = 10 and Ms = mz, 

m| = 4.5 M| - 0.4 M| + 0.3 M2M3 - 0.6 AtM^ -1.3m% + 0.7 + 0.8 mfj^ - 2.0 l/Lip + . . . 

(A.2) 



with that of [78], Eqns. (7) and (11) 

m| = 5.5 iwf - 0.4 + 0.5 M2M3 - 0.8 ItMs - 1-3 + 0.7 + 0.7 fhl^ - 2.2 + . . . 

(A.3) 



or that of [77], Eq. (7) 

m| = 5.2 M| - 0.4 M| + 0.5 M2M3 - 0.8 AjMs - 1-3 fh%^ +0.7 + 0.7 - 2.2 + . . . 

(A.4) 

On the other hand, we agree within 10% with the two-loop result for (1 TeV) of [79], 
Eq. (2.1). 

B Fine-tuning 

To find the least fine-tuned regions of parameter space, we quote again Eq. (2.16), vahd 
for Ms = 1 TeV and tan /? = 20, 

m| = 0.2 Al - 0.7AoMi/2 + 2.9 Mf^^ - 2.1 \i2\^ - 1.3 + 0.7 + 0.8 + . . . 

(B.l) 

A common measure of fine-tuning [80] is derived from the logarithmic sensitivity of the 
electroweak scale with respect to parameter variations, 

0^="^, (B.2) 
a log a 

where a G {Mi/25 ^o, ^/f^' ^Qs' ^f^s' ^^3' • • •} ^^^^ independent dimen- 
sionful GUT-scale parameters. The fine-tuning is then estimated as 

— max Ca . (B.3) 

tine-tumng a 

The worst offender in high-scale mediation scenarios is generically the gaugino mass. This 
remains true even if l^ol is large enough to allow for maximal mixing: From Eq. (2.16) we 
find 

(B.4) 



2.9 


M!/2 


0.35 „ / 




m| 


m| 


0.2 


^0 _ 


0.35 . 




m| 





For l^ol < 3 Ml/25 Cmi/2 dominates. Furthermore, from the fine-tuning point of view, 
it is favorable to go to maximal mixing in order to raise the Higgs mass, rather than to 
raise the overall soft mass scale. For the least fine-tuned regions compatible with a Higgs 
mass m/jO > 123 GeV, with Mi/2/mz ^ 10 and maximal mixing, we find a fine-tuning 
measure of around a few permille. It is clear that the LHC Higgs mass results, when firmly 
established, will substantially raise the fine-tuning price of the MSSM. 
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C Charge- and color-breaking minima 



We briefly review a criterion for charge- and color-breaking minima [81-83] in the context 
which is relevant for us. Consider the direction \Us\ — \Qs\ = \hu\ = X in field space, with 
all other VEVs vanishing. The potential energy along this direction reads 

V{X) = 3 + 2 At yt X^{ml^ + + + X^ + merms . (C.l) 

It is minimized at 




and is negative at Xq if At satisfies 

Al > 3{ml^ + + mj,^ + • (C.3) 

For mz < Ms, the RHS of this inequality is 3{mfj^ + m^^ + mj^^ + ^ 6M|. The 
potential energy of the realistic electroweak vacuum is 

^0 = (g'^ + gl) cos^ 2/3 > -(87 GeV)4 . (C.4) 



If V{Xo) < Vo, there is a charge- and color-breaking vacuum. It is easily checked that, 
for soft terms in the TeV range, the domain V{Xo) < Vq is reached quickly once At starts 
exceeding the bound Eq. (C.3): The electroweak vacuum is quite shallow, compared to the 
CCB vacuum whose scale is set by the soft terms. We therefore conclude that points with 

\Xt/Ms\ > V6 (C.5) 

lead to charge and color breaking. 

Note that this is only a sufficient criterion, and that generally stronger constraints 
can be found by exploring other directions in field space [84]. Furthermore, we are ne- 
glecting i?-terms and loop corrections. Finally, we have evaluated all running quantities 
at the scale Ms here. Demanding the absence of CCB vacua at other scales may, again, 
lead to more restrictive bounds [84]. On the other hand, a CCB vacuum may be viable 
phenomenologically if the lifetime of our false vacuum is long on cosmological timescales. 
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